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Abstract 

Persistence of populations at range edges relies on local population dynamics and fitness, in the case of geographically 
isolated populations of species with low dispersal potential. Focusing on spatial variations in demography helps to predict 
the long-term capability for persistence of populations across the geographical range of species' distribution. The 
demography of two ecological and phylogenetically close macroalgal species with different life history characteristics was 
investigated by using stochastic, stage-based matrix models. Populations of Ascophyllum nodosum and Fucus serratus were 
sampled for up to 4 years at central locations in France and at their southern range limits in Portugal. The stochastic 
population growth rate (X s ) of A. nodosum was lower and more variable in central than in southern sites whilst for F. serratus 
this trend was reversed with X s much lower and more variable in southern than in central populations. Individuals were 
larger in central than in southern populations for both species, which was reflected in the lower transition probabilities of 
individuals to larger size classes and higher probability of shrinkage in the southern populations. In both central and 
southern populations elasticity analysis (proportional sensitivity) of population growth rate showed that fertility elements 
had a small contribution to X s that was more sensitive to changes in matrix transitions corresponding to survival. The 
highest elasticities were found for loop transitions in A. nodosum and for growth to larger size classes in F. serratus. 
Sensitivity analysis showed high selective pressure on individual growth for both species at both locations. The results of 
this study highlight the deterministic role of species-specific life-history traits in population demography across the 
geographical range of species. Additionally, this study demonstrates that individuals' life-transitions differ in vulnerability to 
environmental variability and shows the importance of vegetative compared to reproductive stages for the long-term 
persistence of populations. 
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Introduction 

Understanding the factors that set species geographical borders 
is a central issue in ecology and evolutionary biology [1] that has 
received increased attention in the face of the rapid distributional 
shifts associated with climate change. All species have limits to 
their geographical ranges beyond which they are not found, 
presumably because at these locations their tolerances and 
capacities are constrained or dispersal from source populations is 
limited [2] , but these assumptions have rarely been tested [3] . The 
geographical ranges of species distributions may result from the 
interplay between many factors including habitat availability, local 
environmental conditions, metapopulation dynamics, biological 
interactions and species-specific traits, such as physiological 
tolerances, dispersal potential and vital rates [4], [5]. 



Variations in population structure and dynamics across the 
species distributional ranges must reflect geographical variation in 
the dynamics of vital rates [6], [7]. Therefore, to understand the 
determinants of species ranges it is fundamental to compare the 
vital rates at range boundaries to those of central populations and 
to understand how spatial variation in fitness translates into 
population-level differences in abundance [3]. Such comparisons 
are surprisingly rare, and even less frequent for southern latitude 
margins that are more susceptible to change as a consequence of 
climatic change. Exceptions include empirical comparisons 
between demographic traits at marginal and core locations for 
plants [7], [8], and marine organisms [9], [10]. Studies on 
seaweeds, compared isolated [11] and integrated [12] demo- 
graphic traits at central and southern range locations, or 
developed demographic models for central and marginal popula- 
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tions. However these studies did not include the most marginal 
populations [13]. Understanding the basis of the contrasting 
dynamics of edge and central populations will increase the 
reliability of models on the current, likely threatened e.g. [14], 
[15], status of populations at range margins and to predict their 
long-term response to current and future environmental condi- 
tions. 

Geographically isolated populations of dioecious species with 
short direct dispersal are more likely to be affected by disturbance, 
especially if they are at the margins of their distributional range. 
Among this group of species are the two canopy forming fucoid 
seaweeds Ascophyllum nodosum and Fucus serratus that occur on 
sheltered intertidal shores in the north Atlantic. These are two 
ecologically and phylogenetically close foundation species with 
very different life-spans and phylogeographical histories. A. nodosum 
is a long-lived species able to produce aggregated clonal 
individuals by vegetative sprouting from its holdfast (base), with 
slow growth rate and an estimated maximum life span of a 
genetical individual of hundreds of years [16], whereas F. serratus 
does not produce clonal copies from the holdfast and has a shorter 
life span of about 3 to 4 years [17]. Patterns of genetic diversity at 
the southern range edge are also very different for both species 
[18], [19], [20]. Lower genetic diversity and high between 
population differentiation have been reported for F. serratus 
populations of the Iberian Peninsula when compared with central 
ones [18], [19] which might limit the evolutionary potential of 
such marginal populations [21]. Similar genetic diversity and very 
low differentiation were found across the entire range of 
distribution of A. nodosum in the Adantic, with exception of 
Brittany where higher values were reported [20]. This shows that 
genetic diversity is not reduced nor distinct at the southern edge. 
Population structure and demographic variables of both species at 
marginal and central locations are also different. A. nodosum 
marginal populations exhibited differentiated phenotypic traits in 
relation to central populations with higher values of reproductive 
output and densities and lower individuals' sizes, while F. serratus 
core and border populations did not show differences in vital rates 
[12]. 

The aim of the present study is to compare the demography of 
marginal and central populations of A. nodosum and F. serratus using 
stochastic matrix population models. These models are powerful 
tools to investigate the demography of stage-structured popula- 
tions and allow the incorporation of temporal and spatial sources 
of variability in population dynamics [22]. Environmental 
stochasticity affects population dynamics and evolution of life- 
history traits and thus stochastic matrix models that include 
variations in vital rates due to temporal random fluctuations in the 
environment, were used in this study. Stochastic stage-based 
demographic models for the two species are developed and 
analyzed with focus on stochastic population growth rate, mean 
stable stage distribution, as well as specific sensitivity and elasticity 
of the stochastic population growth rate to changes in stage- 
specific vital rates. More specifically, we hypothesize that the 
stochastic population growth rate should be lower and more 
variable in the marginal than in the central populations and that 
the sensitivity of the stochastic population growth rate should be 
more sensitive to changes in fertility in the marginal populations. 
Furthermore, these analyses aim to contribute to global under- 
standing of demographic contrasts in marginal versus central 
populations and their consequences for persistence of populations 
at their low-latitude distributional limits. 



Materials and Methods 

Study Site 

This study was conducted on the northern coast of Portugal and 
the northern coast of Brittany in France covering, respectively, the 
southern-edge and central populations of F. serratus and A. nodosum. 
These species have maintained stable isolated southern-edge 
populations on the northern coast of Portugal for at least the last 
50 years [23]. In Portugal, the study of F. serratus was performed in 
Viana do Castelo (41°41'27"N; 8°50'57"W) and Amorosa 
(41°38'26"N; 8°49'22"W) and for A. nodosum only in Viana do 
Castelo, the southernmost population in the eastern Atlantic 
(41°41'49"N; 8°51'08"W). In the French coast, two populations of 
both species were sampled, in Roscoff (48°43'38"N; 3°59'22"W) 
and in Santec (48°42'46"N; 4°01'23"W). 

The central populations of both species were located in 
sheltered bays, protected by groups of coastal islands while the 
southern edge populations are on an open coast, exposed to 
prevailing NW swell. However, the area of the shore where A. 
nodosum and F. serratus occur is protected from the direct action of 
waves by a large extension of rocks located in front of their areas. 
This is demonstrated by the prevalence, in these sheltered zones of 
the shore, of brown algae assemblages dominated by Pelvetia 
canaliculata, F. spiralis and F. vesiculosus and A. nodosum. 

No specific permissions were required to perform this study 
since it did not involve collection of specimens and none of the 
species studied had endangered or protected status. 

All sites were strongly dominated by fucoid canopy at mid- 
intertidal sheltered shores. A. nodosum was found in the mid-high 
intertidal zone and F. serratus in the mid-low levels. The maximum 
tidal amplitude is approximately 4.0 m on the northern Portu- 
guese coast and 9.0 m in the northern part of Brittany. Mean 
annual temperatures recorded during the study period under the 
A. nodosum canopy were 1 6°C in Portugal and 1 1 °C in France [12]. 

Sampling Procedure 

For each location (Portugal and France), at each site (Viana do 
Castelo and Amorosa in Portugal, and Roscoff and Santec in 
France) and for each species, three 1 x 1 m quadrats were 
randomly selected and marked in the corners with plastic tags 
and screws. However, for A. nodosum in Portugal only Viana do 
Castelo was sampled since it is the single southern border 
population that exists. The nearest known population occurs 
10 s of kms further north, in Spain. For this reason, twelve 1 x 1 m 
plots were randomly spread along the 2000x30 m total distribu- 
tion of the population in Viana do Castelo. Inside each quadrat 
and at each sampling occasion all the individuals were mapped. 
Individuals with a size large enough to hold a plastic cable tie or a 
glass bead tag were tagged. The basic information needed for the 
size structured matrix model is the fate of individuals, that is, if 
they survive or not. If they survive, then the change in size is 
important in order to calculate the probability of transition of 
individuals from one size class to another. Further, the number of 
recruits is needed to estimate the fertility. For that reason all new 
recruits within each quadrat were registered at each sampling 
occasion. To estimate individual size, the maximum length (1) and 
circumference (c) (in cm) of each individual were measured. The 
predictive equations previously obtained for each species at each 
location [12], for the relationship between algal volume [length of 
the longest frond x maximum circumference squared (lc )] and 
dry weight, were used to calculate the dry biomass (in g) of all the 
individuals measured in the field. Sampling was done from 2005 to 
2009 at variable dates, depending on the knowledge available 
about each species and each location prior to the start of the 
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experiment. For A. nodosum in France, based on previous personal 
observations about the occurrence of a single reproductive period 
yearly, sampling was done once a year from 2006 to 2009 
(summer). Because similar information was not available for F. 
serratus in France, sampling was done twice a year during the first 
year (2006-2007) (winter and summer) showing that the repro- 
ductive period was restricted to few months per year. Thus, 
sampling was done yearly in the following sampling occasions 
(summer). For Portugal, no information about the reproductive 
period of either species was available. A. nodosum'^ receptacles were 
present all year around and thus sampling was performed 3 times 
per year in 2005 (winter, spring and summer). In 2006 sampling 
was done twice a year (winter and summer). Additional in-situ 
observations showed that gamete release occurred mainly between 
October and February, thus between 2007 and 2009 sampling was 
done once a year (summer). For F. serratus in Portugal, sampling 
was done twice a year between 2006 and 2008 (winter and 
summer) and once a year in 2009, since mature receptacles were 
mainly observed between May and September (end of summer). 

For the demographic analyses, individuals were divided into six 
classes: the first class included only the new recruits; individuals in 
this class can only grow to a larger class or die. The other classes 
were based on individual size (<5 g; 5-15 g; 15-54 g; 54-190 g; 
> 190 g), following [24] and where individuals larger than 5 g 
were assumed to be fertile. It was not possible to associate recruits 
to a given fertile alga in the field. For that reason an indirect 
method was used to calculate the fertility values for a given size 
class. The individual fertility values were calculated as contribu- 
tions to the observed number of recruits, for algae in size-classes 3- 
6 in proportion to their biomass. The assumption that fertility is 
proportional to adult biomass is supported by results in earlier 
papers [13], [25], [26], [27]. Projection matrices were constructed 
based on the life cycle graph shown in Fig. 1. 

Model 

The matrix model. Matrix population models are structured 
models where the population is divided into discrete stages, e.g. 
age, size or natural life-cycle stages or a combination of these types 



of stages. The population matrix describes the transition proba- 
bility from one stage to another and the probability to stay in the 
same stage as well as the stage specific fertility. The change in 
population size can be projected in discrete time steps and there 
are numerous analyses that can be applied to the model (see [22] 
for a review). A stochastic matrix model is usually based on a finite 
set of projection matrices {Ao, A b A k }, each representing a 
specific environment (e.g. here year and country). The stochastic 
process can be described by a first-order Markov chain of length T 
(t=0, 1, Tj, where Tis a large number (usually thousands of 
time steps). To model the population growth a matrix A, is drawn 
from the subset of matrices at each time t, with a set probability 
and multiplied by a column vector a(t) specifying the number of 
individuals in each stage class at time t. The basic stochastic matrix 
population model can be expressed as 



n(r+l) = A ( n(0 



(1) 



where n(<) is a population vector and A t is the population 
projection matrix generated by the stochastic environmental 
process at time step t. If the probability distribution of the 
different environments does not change through time the 
environment is considered homogenous. The stochastic popula- 
tion growth rate can then be calculated as 



logA, 



lim — log||A r _ 



■ A 0 n 0 || 



(2) 



To our knowledge there are no statistical methods to compare 
the population growth rate and thus is not possible to calculate the 
uncertainty in the estimation of this parameter. 

Model analyses and simulations. Due to the unbalanced 
temporal sampling there are three matrices describing the first 
year and two for the second year for A. nodosum in Portugal. In year 
three and four there is only one matrix per year similarly to A. 




\ Mortality (M) \^ Probability of staying in the same size class (L) 
t Fertility (F) | Transition probability to bigger size classes (G) 



; Transition probability to smaller size classes (S) 



Figure 1. Life cycle graph for A. nodosum and F. serratus showing the different stage classes and all possible transitions. Individuals in 
a size class (2-6) that survive can grow (G), shrink (S) or stay in the same size class (L). Class 1 includes new recruits that, if survive, enter class 2. 
Individuals in class 3-6 are reproductive. 
doi:10.1371/journal.pone.0092177.g001 
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nodosum in France where there is only one matrix for each of the 
three years. For F. serratus in Portugal there are two matrices per 
year for the first and second year and one matrix for the third year. 
For F. serratus in France there are two matrices for the first year 
and one matrix per year for years two and three. For years with 
more than one matrix a periodic model was applied to calculate a 
yearly matrix [28]. The stochastic population growth rate was 
calculated using yearly matrices where all years had the same 
probability to appear in the stochastic sample path. For each 
species and site a mean matrix was calculated from the yearly 
matrices (Table 1) e.g. for a site with three yearly matrices the 
mean matrix A m = (A!+A 2 +A 3 )/3. The stable stage distribution 
and reproductive values were calculated as the right and left 
eigenvectors corresponding to the dominant eigenvalue of the 
mean matrix for each species and country. 

Sensitivity and elasticity. The sensitivity and elasticity 
(proportional sensitivity) of the stochastic population growth rate 
to changes in the matrix elements were calculated using model 
simulation. Due to the unbalanced temporal sampling the methods 
in [29] and [30] could not be used directly but the simulations are 
equivalent to these methods. To calculate the sensitivity of the 
stochastic population growth rate a small figure (0.0001) was 
added to one element at the time in all the basic matrices. For 
years with more than one matrix, the yearly matrix was then 
calculated with the periodic model and added to the set of yearly 
matrices. A new stochastic population growth rate was calculated 
for the perturbed set of matrices using the same stochastic sample 
path as when calculating the original growth rate. This new 
stochastic population growth rate minus the original one divided 
by the size of the perturbation -1 gave the sensitivity for the matrix 
element that was perturbed. The matrices were then set back to 
their original values and the next element was perturbed until 
values for all elements were obtained. The size of the perturbation 
was decreased until the sensitivity values converged to three 
significant digits. The elasticity of the population growth rate was 
calculated in a similar way but with a factor close to unity (1.0001) 
instead of an absolute perturbation. All calculations and simula- 
tions were performed in MATLAB R2011a. 

Results 

Stochastic Population Growth Rate and Stable Stage 
Distribution 

Model simulations showed that the stochastic population growth 
rate (X^j varied between species and locations. When X s is > 1 the 
average population size is growing and when it is < 1 the average 
population size is decreasing. X s was higher for A. nodosum in 
Portugal (0.97) than in France (0.87) but lower for F. serratus in 
Portugal (0.77) than in France (1.05) (Fig. 2). All the X s were 
slightly below 1, except for F. serratus in France (Fig. 2). For A. 
nodosum, variation in population growth rate was larger in France 
than in Portugal but the opposite trend was recorded for F. serratus 
(Fig. 2). 

The stable stage distribution calculated from the mean matrix 
showed for both A. nodosum and F. serratus, that populations in 
France had higher proportion of individuals belonging to the 
largest size classes (classes 5 and 6 for A. nodosum and classes 3 to 6 
for F. serratus) than in Portugal. However, all the populations, 
especially F. serratus, both in Portugal and France, were dominated 
by smaller individuals belonging to classes 1 and 2 (Fig. 3). 

Reproductive Values 

Reproductive values increased with plant size for both species at 
both locations but more pronouncedly for F. serratus, especially in 



Portugal (Fig. 4). These results were related to the fertilities of each 
class and the transition probabilities between size classes. Higher 
fertilities were found in Portugal for both species at all size classes. 
However, for A. nodosum in Portugal the probability of transition to 
the class with higher fertility (class 6) is very low and the shrinkage 
to smaller size classes much higher than in France (Table 1). On 
the contrary, for F. serratus the differences in the transition to 
higher size classes and in the shrinkage of individuals to smaller 
sizes between locations are not so pronounced, and do not balance 
the much higher values of fertility found in Portugal (Table 1). 

Elasticity and Sensitivity 

The pattern of elasticities of the stochastic population growth 
rate to changes in the matrix elements showed that the loop 
transition was the most important for the X s of A nodosum while the 
growth transition was the most important for the X s of F. serratus 
(Fig. 5). Elasticities were very similar for F. serratus populations in 
France and in Portugal, with the shrinkage transition as the least 
important for X s (Fig. 5). Fertility transitions were very low for A. 
nodosum when compared to F. serratus but were slighdy higher in 
Portugal than in France (Fig. 5). For A. nodosum, loop transitions 
were more important in Portugal than in France while growth and 
shrinkage transitions showed the opposite geographical pattern 
(Fig. 5). The comparison between the elasticity across size classes, 
without taking into account the fertility, showed that for A. nodosum 
in France there was an important contribution of the 3 largest size 
classes (4—6) to X s while in Portugal it was classes 2 to 4 that 
influence X s the most (Fig. 6). For F. serratus the contribution of 
each size class to X s was more evenly distributed among size classes 
1 to 5 for both locations (Fig. 6). 

The pattern of sensitivities of the stochastic population growth 
rate to changes in the matrix elements showed that, globally, there 
was a selection for individual growth in both species at both 
locations. However, for A. nodosum the sensitivity to changes in loop 
and shrinkage transitions were also important for population 
growth, while the sensitivity to changes in growth transition was 
totally dominant over the other matrix elements in F. serratus 
(Fig. 7). In general, the sensitivity of the stochastic population 
growth rate to changes in fertility was very low, in particular for F. 
serratus. For A. nodosum, fertility transition has higher sensitivity in 
Portugal than in France (Fig. 7). For A. nodosum the sensitivity to 
changes in growth and loop transitions were higher in Portugal 
while the shrinkage transition was slighdy more important in 
France. The comparison between the sensitivities, without taking 
into account the fertility across size classes, showed that the 
sensitivity of the stochastic population growth rate was higher for 
changes in classes 2 to 4 for A. nodosum in Portugal, while in France 
it was in general lower and balanced between classes 2 and 6 
(Fig. 8). For F. serratus sensitivity of X s was much higher for changes 
in smaller classes (1 and 2) and higher in Portugal than in France 
(Fig. 8). 

Discussion 

This study revealed distinct long-term demographic trends for 
central and edge populations of A. nodosum and F. serratus, 
demonstrating the role of species identity and related life-history 
traits on the dynamics of edge populations. These results have 
wide implications for the understanding and predicting of 
biodiversity changes at range edges, by revealing that marginal 
versus central populations can have important demographic 
differences that determine their capacity for persistence or 
likelihood to shift their distributional limits. 
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Figure 2. Mean (±SD), minimum (*) and maximum (#) 
stochastic growth rates (>. s ) for southern (Portugal) and central 
(France) A. nodosum and F. serratus populations. 
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A. nodosum and F. serratus are contrasting species in terms of life 
histories, habitats and genetic background. Earlier studies show 
that phenotypic traits are differentiated between A. nodosum'':, 
southern-edge and central populations while F. serratus has similar 
traits at both locations but occupies a narrower vertical 
distribution in the south [12], [21]. The results of the present 
study show how these differences in fitness related traits are 
translated into the demography of the species and how sensitive 
the stochastic population growth rate is to changes in vital rates. 
An extension of this could be to use these models to predict 
population persistence under the future scenario of environmental 
variability posed by climate change. 

Stochastic Population Growth Rate and Stable Stage 
Distribution 

For A. nodosum, modeled X s was lower and more variable at the 
center than at the edge of distribution, but below 1 for both 
locations. For F. serratus X s was much lower and variable at 
southern than at central locations. In general, macroalgae 
populations have lower and more variable X s in environmentally 
variable habitats as a consequence of high mortality of individuals 
[26], [31]. An increase in the mortality of A. nodosum individuals 
has been related to different sources of disturbance such as ice 
scouring, harvesting and trampling [26], [32], [33]. The lower X s 
found for the A. nodosum population in France was mainly the result 
of the lower survival of individuals in all size classes when 
compared to Portugal (see Table 1). During this study, deposition 
of sand was frequently observed in some of the study plots in 
France but never in Portugal. This episodic disturbance might 
have accounted for the lower survival rates observed in France. A 
previous study in the same area of France in the late 1990's found 
much lower values of mortality and never reported accumulation 
of sand in the study plots, but still found X s of less than 1 due to low 
recruitment rates (Aberg et al. personal observation). Values of X s 
below 1 were reported for A. nodosum in the northern part of the 
Swedish west coast, close to its distributional limit further south on 
the west coast [13]. For central populations, X s varied from slighdy 
below to slightly above 1 , depending on the location and sampling 
dates [13], (Aberg et al. personal observation). This demonstrates 
the need to conduct long-term observational studies covering a 
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comprehensive spatial extent, along the distributional range of 
species, when constructing long-term predictions about popula- 
tion's trends. 

Southern range F. serratus 's populations show reduced plasticity 
to environmental variation [21]. In the present study, model 
simulations showed, for F. serratus, that environmental fluctuation 
at the boundary of its distribution, might cause its disappearance 
from previously colonized areas, which take a long time to recover 
[34]. The X s values found in this study might reflect the 
maladaptation of populations to years with high environmental 
variation. 

For both species, large individuals were much more abundant in 
France than in Portugal. The individual growth rate of both 
species is similar at both locations (Araujo et al. personal 
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Figure 5. Elasticity of the population growth rate to changes in matrix elements summed over different regions (growth, loop, 
shrinkage and fertility) for southern (Portugal) and central (France) A. nodosum and F. serratus populations. 

doi:1 0.1 371 /journal.pone.00921 77.g005 

observation), explaining that the lower transition probability to 
larger size classes in Portugal is caused by higher breakage of 
fronds. In previous studies, increased breakage of macroalgae's 
fronds has been related with several factors such as grazing [35], 
[36], [37], [38], wave action [39], [40], harvesting [41] and 
trampling [33]. However, the comparative quantitative effects of 
most of these factors on breakage of the fronds at central and 
southern locations are unknown for these species. 

Reproductive Values 

The reproductive values increased with individual size for both 
species but more pronouncedly for F. serratus, especially in 
Portugal. This indicates that F. serratus individuals in Portugal 
have a higher fitness advantage in terms of reproduction and 
survival if they grow to larger size classes. An increase in the 
reproductive biomass with size of the plant is a common trend for 
seaweeds [25], [43] and an increase in the reproductive value with 
size class has been reported in previous studies [24], [31], [41]. 
Reproductive effort might increase following disturbance [44], 
[45] , [46] . The high reproductive values of F. serratus in this study 
might reflect the levels of environmental disturbance found at its 
edge of distribution. 
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Figure 6. Elasticity of the population growth rate to changes in 
the survival transitions for each stage class in southern 
(Portugal) and central (France) A. nodosum and F. serratus 
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Figure 7. Sensitivity of the population growth rate to changes in matrix elements expressed as relative contributions and summed 
over different regions (growth, loop, shrinkage and fertility) in southern (Portugal) and central (France) A. nodosum and F. serratus 
populations. 
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Elasticity 

Elasticity of the population growth rate to changes in survival 
was higher than for changes in fertility, for both species at both 
locations. A. nodosum was more susceptible to changes in loop 
transition followed by growth in both areas, while changes in 
growth transitions were the most important for F. serratus followed 
by loop. However, F. serratus populations were more susceptible to 
changes in fertility rates than those of A. nodosum. In general for 
algae, higher elasticity values for survival than for fertility are 
reported for longer-lived species, while in shorter-lived species 
elasticity values for fertility, although lower than for survival, are 
generally higher [13], [26], [32], [48], [47]. 

The higher susceptibility to changes in loop transitions of A. 
nodosum populations demonstrates that their long-term persistence 
might be affected by factors increasing the shrinkage of individuals 
(in particular of adults). Among these factors, an increase in the 
human frequentation of coastal areas might be important since 
human trampling negatively affects the communities at the 
southern-limit of A. nodosum, both in the short and long term 
[33]. The density of invertebrates inhabiting A. nodosum adults is 
lower in Portugal than in central and northern populations [12], 



[36]. However, in spite of the decreasing trend in limpet grazing 
from north to south [42], the opposite trend was recorded under 
the A. nodosum canopy when comparing edge and central 
populations [12]. Harvesting, although not important in Portugal, 
occurs in France, and increases individual shrinkage [32], [49]. 
Also changes in the frequency of storms predicted by models of 
climate change [50] or in the pressure of grazing [51] might 
account for increased shrinkage of individuals. Although the A. 
nodosum populations in Portugal were slighdy more susceptible to 
changes in fertility than in France, elasticities for fertility 
transitions were low. The reproductive output [12] and reproduc- 
tive allocation (Araujo et al. personal observation) of A. nodosum 
was higher in Portugal than in France, despite the apparent low 
importance for population growth and persistence demonstrated 
by the model. However, in long-lived algae, recruitment is 
important for the long-term stability of populations [41], [52], 
and massive recruitment may occur under exceptional circum- 
stances such as disturbance events [31], [53]. Short-term studies, 
like the one presented here, are likely to miss these extraordinary 
recruitment events, potentially contributing to the occurrence of 
estimates of population growth below 1 observed in this study. 
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Figure 8. Sensitivity of the population growth rate to changes 
in the survival transitions for each stage class in southern 
(Portugal) and central (France) A. nodosum and F. serratus 
populations. 
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Factors affecting growth of individuals also have negative effects 
on the growth of the population, as demonstrated by the high 
elasticities of growth transition found for both species. The 
documented North Atlantic trend of raising sea and air 
temperatures [54] is expected to have more immediate and severe 
consequences at species borders. As expected, temperatures in the 
intertidal area are higher in Portugal than in France [12]. 
Increased sea surface temperature will decrease A. nodosum growth 
at temperatures above 19°C with lethal effects shown above 28°C 
[55]. F. serratus southern edge populations show much lower 
resilience to heat stress than central ones [21]. Consequendy, it 
can be expected that, under the current scenario of global 
warming, these southern limit populations will suffer a decrease in 
individual growth that could have important consequences for 
population persistence. 

Sensitivity 

The sensitivity of the stochastic population growth rate to 
changes in a given vital rate can be interpreted as the selection 
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